Search for muon-neutrino to electron-neutrino transitions in MINOS 
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This letter reports on a search for i^^ transitions by the MINOS experiment based on a 

3.14 X 10^" protons-on-target exposure in the Fermilab NuMI beam. We observe 35 events in the 
Far Detector with a background of 27 ± 5(stat.) ± 2(syst.) events predicted by the measurements 
in the Near Detector. If interpreted in terms of — > oscillations, this 1.5 a excess of events 
is consistent with sin (26I13) comparable to the CHOOZ limit when |Am^|=2.43x 10"^ eV^ and 
sin^ (2^23) = !. are assumed. 

PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 



Several experiments have provided compelling evi- 
dence for muon neutrino disappearance as a function of 
neutrino energy and distance traveled [H, 0, H, 0, H, @, 01 ■ 
These observations support the description of neutrinos 
in two distinct mass and flavor bases, related by the 
3x3 neutrino mixing matrix Q. The MINOS experi- 
ment provides the most precise measurement of the atmo- 
spheric mass splitting, IAto^] = (2.43±0.13) x 10"^ eV^ 
[3, At this mass scale, the dominant oscillation 

channel is expected to be Vy^ Vt, but sub-dominant 
j/^ — > Ve transitions are not excluded [lo| . Observa- 
tion of appearance would imply a non-zero value of 
6'i3, opening the possibility of observing CP violation 
in the leptonic sector. The current best experimental 
limit Ji.], implies sin^ (26*13) < 0.15 at the 90% confi- 
dence level (C.L.) for the MINOS jAm^j value. In ad- 
dition to these parameters, the probability of electron 
neutrino transitions in MINOS depends on sin^(26'23), 
the CP violation parameter, 5c p, and the sign of Am?. 
Two other experiments have given limits with less sen- 
sitivity [HillSl; MINOS is the first experiment to probe 
sin'^(20i3) with sensitivity comparable to the CHOOZ 
limit at |Am2 1=2.43 X 10-3 and sin^ (26'23)=1.0. 

Neutrino interactions from a beam produced by the 
Fermilab NuMI facility [11] are recorded at the MINOS 
Near (ND) and Far (FD) Detectors, located at 1 km and 
735 km, respectively, from the production target. High 
statistics data from the ND establish the properties of 
the beam before oscillations. Observation of additional 
j/e interactions in the FD relative to the ND provides ev- 
idence of Vy //g oscillation. The two detectors are of 
similar design to reduce systematic uncertainties from the 
physics of neutrino interactions, the neutrino flux, and 
detector response 0, [iM [3 ■ The detectors are magne- 
tized, tracking calorimeters composed of planes each with 
layers of 2.54 cm thick steel and 1.0 cm thick scintillator 
(1.4 radiation lengths per plane). The scintillator layer 
is composed of 4.1 cm wide strips (1.1 Moliere radii). 

The beam is comprised of 98.7% v^+Vfj_ and 1.3% 
t'e+i^e- The latter originate from decays of muons pro- 
duced in pion decays and from kaon decays. The v^. flux 
below 8 GeV is largely from muon decay and is well con- 
strained by the measured flux 0, [l3l- The present 
analysis is based on an integrated exposure of 3.14 x 10^" 
protons delivered to the NuMI target. 



The search for electron neutrino appearance relies on 
identifying charged current (CC) i^e + Fe e+X interac- 
tions that produce an energetic electron [I^. This elec- 
tron initiates an electromagnetic cascade and deposits its 
energy in a relatively narrow and short region in the MI- 
NOS calorimeter. Additional calorimeter activity is pro- 
duced by the breakup of the recoil nucleus (X). Other 
neutrino scattering processes can produce similar event 
topologies in the MINOS detector. These include neutral 
current (NC) v + Fe + X interactions and t'^-CC in- 
teractions with low energy muons, both having hadronic 
showers with an electromagnetic component arising from 
7r° decays. Less significant backgrounds arise from in- 
trinsic beam i^e-CC interactions, Vt-CC interactions from 
oscillations, and cosmogenic backgrounds. 

We select events with reconstructed energy between 1 
and 8 GeV, encompassing the maximum of the — > Vg 
oscillation probability. The lower limit mainly removes 
NC events, while the higher limit removes beam i/g-CC 
events from kaon decays. Additionally, events are se- 
lected to be in time with the accelerator beam pulse, 
and directional requirements are applied to limit back- 
ground from cosmogenic sources to less than 0.5 events 
(90% C.L.). Events are required to have a reconstructed 
shower and at least 5 contiguous planes each with energy 
deposition above 1 MeV. Events with tracks longer than 
25 planes are rejected. Monte Carlo (MC) simulations in- 
dicate that these cuts improve the signal-to-background 
ratio from 1:55 to 1:12, assuming sin^ (26'i3)= 0.15. 

Further enrichment of the i/g-CC selected sample is 
achieved using a method based on an artificial neural 
network (ANN) with 11 input variables characterizing 
the longitudinal and transverse energy deposition in the 
calorimeter that separate the signal i/g-CC events from 
NC and J/^-CC background [l^ ■ The acceptance thresh- 
old is determined by maximizing the ratio of the accepted 
signal to the expected statistical and systematic uncer- 
tainty of the background. With these selection criteria, 
and assuming sin^ (26*13)= 0.15, this method gives a 1:4 
signal-to-background ratio. 

The Monte Carlo simulation of the beam line and de- 
tectors is based on GEANT3 [2^ and the hadron produc- 
tion yields from the target are calculated by FLUKA 21 1. 
The calculated neutrino flux is adjusted to agree with the 



ND Vfj^-CC data [17|. Neutrino interactions and further 
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FIG. 1: Reconstructed ND energy spectra of the i/e-CC se- 
lected backgrounds from NC (dashed) and t'^-CC (dotted) in- 
teractions as obtained from the horn-off method. The shaded 
histogram shows the beam Ve-CC component from the simu- 
lation. The solid histogram corresponds to the total of these 
three components which are constrained to agree with the 
data points. The uncertainties shown on the data are sta- 
tistical and are not visible on this scale; uncertainties on the 
components are systematic. 



re-interactions of the resulting hadrons within the nu- 
cleus are simulated with NEUGEN3 [13] ■ Uncertainties in 
the composition and kinematic distribution of the parti- 
cles that emerge from the nucleus can be large, but these 
and many other uncertainties mostly cancel when com- 
paring neutrino interactions between the ND and FD. 

Events selected by the ANN in the ND predict the 
number of background events expected in the FD and 
reduce reliance on the simulation. Considerations of 
oscillations and beam line geometry require that each 
background component, v^-CC, NC, and beam i/g-CC, 
be treated separately in the prediction of the FD back- 
grounds. The first two background components are de- 
termined using an NC-enriched data sample recorded 
with the NuMI magnetic horns turned off. In this config- 
uration the pions are not focused; the low-energy peak of 
the neutrino energy distribution disappears, leaving an 
event sample dominated by NC events from higher en- 
ergy neutrino interactions. These data are used in con- 
junction with the standard beam configuration data, and 
the simulated ratios of the horn-on to horn-off rates for 
each component, to extract the individual NC and i/^- 
CC background spectra. The smaller beam z^e-CC com- 
ponent is calculated from the MC using fc'p-CC events 
observed in the ND. 

Figured] shows the data in the ND and the derived NC, 
v^-CC, and beam z/g-CC backgrounds. The ND back- 
ground is (57±5)% NC, (32±7)% z/^-CC and (11±3)% 



beam Ve-CC events. Systematic errors on the compo- 
nents arise from uncertainties in the beam flux, cross- 
section and selection efficiency. The errors on the NC 
and i^^-CC components are derived from the data and 
are correlated due to the constraint that the background 
must add to the observed ND event rate. The uncer- 
tainty on the beam i^e-CC comes from the z^^-CC events 
observed in the ND . 

As a crosscheck, a second technique to study the ND 
background sample uses an independent sample of show- 
ers from f^-CC events selected with long tracks [7|. The 
hits associated with the muon track are removed from 
the event, and the remnant showers are subsequently an- 
alyzed as a sample of NC-like events [1^ [IS] • The pro- 
cedure is applied to the data and MC, and the i/g-CC 
selections are applied to both. Differences between the 
muon-removed data and muon-removed MC samples are 
used to adjust the predicted NC background. As in the 
first method, the beam i^e-CC background is taken from 
the MC, and the remainder of the observed ND back- 
ground are classified as v^-CC events. The background 
components calculated from the muon-removed sample 
agree with those obtained from the horn-off method [25| . 

After decomposition of the ND energy spectrum into 
background components, each of these spectra is multi- 
plied by the Far to Near energy spectrum ratio from the 
simulation for that component, providing a prediction of 
the FD spectrum. The simulation takes into account dif- 
ferences in the spectrum of events at the ND and FD due 
to beam line geometry as well as possible differences in 
detector calibrations and event shape. Oscillations are 
included when predicting the v^-CC component. We ex- 
pect 26.6 background events, of which 18.2 are NC, 5.1 
are v^-CC, 2.2 are beam z/g and 1.1 are Vr [26| . 

To estimate the efficiency for selecting Ve-CC events, 
we use the muon-removed events from data and MC, then 
embed a simulated electron of the same momentum as 



the removed muon. Test beam measurements 27| indi- 
cate that the selection efficiency of single electrons agrees 
with the simulation to within 2.6%. The v^-CC selection 
efficiency obtained from the data agrees with the selec- 
tion efficiency obtained from the MC to within 0.3%. We 
estimate our efficiency for selecting Vg-GG events to be 
(41.4±1.5)% [11]. 

Systematic uncertainties are evaluated by generating 
modified MC samples and by quantifying the change in 
the number of predicted background events in the FD. 
Table H] shows that the dominant uncertainties arise from 
Far/Near differences: relative energy scale calibration 
differences (a) , details of the modeling of the photomulti- 
plier tube (PMT) gains (b) and crosstalk (c), and relative 
event rate normalization (d). Other uncertainties result- 
ing from neutrino interaction physics, shower hadroniza- 
tion, intranuclear rescattering, and absolute energy scale 
uncertainties (e) affect the events in both detectors in 
a similar manner and mostly cancel in the extrapola- 
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background events 


Far/Near ratio: 




6.4% 


(a) Relative Energy Scale 


3.1% 




(b) PMT Gains 


2.7% 




(c) PMT Crosstalk 


2.2% 




(d) Relative Event Rate 


2.4% 




(e) All Others 


3.7% 




Horn-off (systematic) 




2.7% 


Horn-ofF (statistical) 




2.3% 


Total Systematic Uncertainty 


7.3% 



TABLE I: Systematic uncertainty in the total number of back- 
ground events in the Far Detector. 



tion. The individual systematic uncertainties are added 
in quadrature along with the uncertainty arising from the 
background decomposition in the ND to give an overall 
systematic uncertainty of 7.3% on the expected number 
of background events selected in the FD. The expected 
statistical uncertainty is 19%. 

The prediction and uncertainty of the backgrounds 
in the FD are established before examining the FD 
data. Additionally, the independent and signal-free 
muon-removed FD data sample is examined. In that 
sample, we observe 39 events, with an expectation of 
29 ± 5(stat.) ± 2(syst.). The selected events were investi- 
gated and no evidence of abnormalities was found. 

Figure [5] shows the FD data as a function of the ANN 
selection variable. The signal acceptance threshold was 
optimized prior to examination of the FD data to be 
0.7. We observe 35 events in the signal region with a 
background expectation of 27±5(stat.)±2(syst.). In the 
region of the selection variable well below the acceptance 
threshold (< 0.55), we observe 146 events, compared to a 
pure background expectation of 132± 12(stat.) ±8(syst.). 
The observed energy spectrum for the events in the signal 
region is shown in Figure [3] 

A second selection method. Library Event Matching, 
is used as a crosscheck. In this technique, each candidate 
is compared to a large library of simulated Ve-CC and 
NC events [2^. This method gives a better background 
rejection than the ANN algorithm, but with increased 
sensitivity to some uncertainties. As in the ANN method, 
we observe a small excess (< 2(7) for the muon-removed 
sample and in the region below the selection cut. In 
the signal region, we observe 28 selected events, with a 
background expectation of 22 ± 5(stat.) ± 3(syst.); these 
results are consistent with the ANN selection. 

Figure m shows the values of sin^ (26*13) and Sep that 
give an excess of events consistent with our observa- 
tion from the ANN selection. The oscillation proba- 
bility is computed using a full 3-flavor neutrino mix- 
ing framework that includes matter effects [2^, which 
introduces a dependence on the neutrino mass hierar- 
chy (the sign of Am^). The MINOS best fit values of 
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FIG. 2: Distribution of the ANN selection variable for events 
in the FD. Black points show data with statistical error bars. 
The non-shaded histogram shows the background expecta- 
tion. The shaded region shows the additional i/e-CC events 
required to explain the observed excess with the oscillation 
hypothesis. 
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FIG. 3: Reconstructed energy distribution of the i^e-CC se- 
lected events in the FD, with the exception of the energy cut. 
No data events are found below 1 GeV or above 8 GeV, consis- 
tent with expectation. Black points show data with statistical 
error bars. The non-shaded histogram shows the background 
expectation. The shaded region shows the additional Ve-CC 
events required to explain the observed excess with the oscil- 
lation hypothesis. 



IAto^ 1=2.43 X 10^3 eV^ and sin^ (26*23) =1.0 are used as 
constants in the calculation fsoj. Statistical fluctuations 
(Poisson) and systematic effects (Gaussian) are incorpo- 
rated via the Feldman-Cousins approach [31] which de- 
termines the confidence intervals. 

In conclusion, we report the first results of a search 
for I'e appearance by the MINOS experiment. The 
35 events in the Far Detector after Ve selection for 
3.14 X lO'^" protons-on-target are 1.5 cr higher than the 
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sin'(2ei3) 

FIG. 4: Values of sin'^ (2^13) and 6cp that produce a num- 
ber of events consistent with the observation for the normal 
hierarchy (top) and inverted hierarchy (bottom). Black lines 
show those values that best represent our data. Red (blue) 
regions show the 90% (68%) C.L. intervals. 

background expectation of 27±5(stat.)±2(syst.). Assum- 
ing |Am2|=2.43xl0-3eV2 and sin^ (26'23) = 1.0, the best 
fit for the normal hierarchy is just below the CHOOZ [llj 
limit for all values oi Sep- 
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